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We discuss the general mechanism leading to long-range rapidity and angular correlations 
produced in high energy collisions (the "ridge"). This effect naturally appears in the high 
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counting and stress the relevance of Pomeron loops. 
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1. Introduction 

While built as a discovery machine, the LHC at CERN has a very rich research 
program exploring QCD in extreme conditions. This program ranges from the fun- 
damental measurement of total scattering cross section in proton-proton (pp) colli- 
sions at very high energies to production and study of Quark Gluon Plasma (QGP) 
in heavy ion collisions. 

One of the first and intriguing observations made by the CMS collaboration^ is 
a discovery of long-range two-particle correlations in pp collisions at y/s = 7 TeV. 
The correlations are observed between two charged hadrons which emerge from 
the collision collimated in their azimuthal angle 4> while being widely separated in 
(pseudo-) rapidity 77. Remarkably, while not present in the minimal bias events, 
the correlations are observed in high multiplicity events. These are extremely rare 
events (approximately 10~ 6 of the total event sample) in which the total multiplicity 
of produced charged hadrons is many (7-8) times larger than a typical multiplicity 
in minimal bias collisions. 

Correlations, similar to the ones observed by the CMS, were seen before at 

191*3141 

RHIC in collisions of heavy ion a 1 1 I The phenomenon was dubbed "ridge" due to 
its shape on the 4>—rj data plot (see below). An accepted explanation of the "ridge" 
at RHIC is related to collective flow in QGP. The fact that a similar ridge structure 
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has been observed in pp collisions may hint at an intriguing possibility that high 
multiplicity pp events produce a state somewhat similar to QGP. This question is 
however still open today, as much more data would be needed to provide a clear 
answer. 

Recently, CMS released its preliminary data from the test run of proton- lead col- 
lisions at 5.02 TeV per nucleorP. The "ridge"-like structure of di-hadron azimuthal 
correlations is again clearly seen in high multiplicity events, and its magnitude is 
significantly stronger than in the pp data. 

There is a theoretical basis for expectation that at asymptotically high energies 
hadrons and nuclei behave alike. At high energies both protons and heavy ions 
appear as dense clouds of gluons. The growth of gluon densities with energy is 
governed by the QCD evolution equations, which are universal for all hadrons. Thus, 
from this point of view, pp collisions at very high energies are expected to resemble 
heavy ion collisions at somewhat smaller energy. It is however overly optimistic to 
expect that this regime has already been reached at LHC, and it certainly has not 
been reached for the minimum bias pp collisions. The rare high multiplicity events 
may represent the tail of the distribution where high density has been achieved 
already, but this too is far from certain. 

The long range rapidity correlations may not be (and probably are not) related 
to the formation of QGP. However they are a striking feature of the observed data, 
and are very likely a universal property of a highly energetic QCD collision, pro- 
vided the number of particles produced in the final state is large enough to probe 

such correlations. Since the correlations extend over a very large rapidity inter- 
nal 

val, causality considerations- (similar to cosmological arguments behind baryonic 
acoustics and CMB fluctuations) lead one to conclusion that the origin of the cor- 
relations is at very early stage of the collision, and in fact most likely in the initial 
state. Without such early time correlations, relativistic particles separated by many 
units of rapidity flying out of the interaction region would be causally disconnected 
and uncorrelated. At the same time, high multiplicities of produced particles must 
be due to high gluon densities of colliding hadrons. We are thus naturally lead to 
a conclusion that the observed correlations are rooted in the high density QCD 
(hdQCD). The nature of this relation is the subject of this review paper. 

There have been numerous publications pr oposing various m echanisms respon- 
sible for the "ridge" in AA and in pp coilision d^ l 9 | 10 | ll | 12 | 13 | It is not our goal 
to provide a full coverage of those. Here we focus on correlations that originate in 
the initial conditions, or the wave-functions of colliding hadrons and their relation 
to high density effects. The main ideas reviewed here, appeared first in Refs.lfMj 
[TBI IT6| . We will explain how the theory of high density QCD naturally leads to 
the correlations of the kind observed in the pp experimental data. The mechanism 
discussed here does not refer at all to collective flow in the final state, and thus 
is not likely to be relevant for the explanation of the ridge in heavy ion collisions. 
We take in this review the point of view that although the origin of long range ra- 
pidity correlations is similar in heavy ion and proton-proton collisions, the angular 
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collimation in the two systems is probably of a different origin. 

IidQCD manifests itself in the phenomenon of gluon density saturation!!^! Gluon 
densities grow with energy. The lo w (Bjor ken) x evolution of gluon densities is gov- 
erned by the linear BFKL equatiorPSHD which leads to a rapid (exponential with 
rapidity) growth. Once the densities reach the "critical" value, where the density of 
gluons per unit phase space is of order of the inverse of the strong coupling constant, 
non-linear effects in the evolution become important. These effects slow down the 
growth of the density, leading to gluon saturation, characterized by the appearance 
of a new scale Q s (x), known as the saturation momentum. Such a dense glu onic 
state is frequently referred to as the Color Glass Condensate 

(CGcjsnElEa The 

strong coupling constant a s , evaluated at the saturation scale Q 2 ~ 5 GeV 2 for 
the LHC energies, is rather small, which makes many aspects of saturation pertur- 
batively accessible. However, even for a small coupling constant, the fact that the 
gluon densities are high invalidates naive perturbation theory. The proper treatment 
of the saturation regime therefore requires resummation of the naive perturbation 
theory. 

The nonlinear general ization of th e BFKL e quation that resums the hi gh density 
effects is known as the Bk P^M JIMWIjI jZ51ZZlZ ai9 | 30 | 31 | 20 | 21 | 22 |^ ctioiial 

equation. This (functional) equation is strictly speaking applicable when one of the 
colliding particles is dilute, such as in the deep inelastic scattering (DIS) explored 
at HERA, or in proton-nucleus collisions. Although this is not the case when both 
colliding particles are high energy hadrons, and the QCD description of such systems 
is still not well understood, we will be using the BK-JIMWLK approach as a proxy 
to qualitative description of hadronic high energy scattering. 

Some aspects of relation of saturation physics with the ridge phenomenon have 
been reviewed in Refs.flMJ 131 GHE [3S 133) ■ 

2. Experimental data on di-hadron correlations 
2.1. The "Ridge" in heavy ion collisions. 

The ridge was first observed in gold-gold collisions at y/s — 200 GeV at RHIC by 
the STAR experiment!^ (Fig. [IJ. The initial measurement of di-hadron correlations 
was performed with a high momentum trigger particle pr > 4 GeV, but subse- 
quently the correlations were measured for both, trigger and associated momenta 
in the (semi) soft range. Here <fi is the azimuthal angle, A</> is the angular difference 
between the momentum of the measured hadron and the direction of the trigger. 
The peak around A(j) — corresponds to the trigger jet fragmentation region. The 
bump at Atfi = n is the usual back-to-back correlation. r\ ~ ln(p z /p±) stands for 
(pseudo-)rapidity of produced particles. For two hadrons with approximately equal 
transverse momenta, A77 measures their relative momenta along the beam direction, 
or, in other words, relative energies. Centrality dependence (indicated in percents 
above each plot) of the "ridge" clearly demonstrates the ridge in central collisions 
while it disappears towards peripheral region. Further studies by PHOBOsP and 
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STAR,El(Fig. |5J) were performed with and without hard trigger and confirmed the 
appearance of the ridge and its essential (semi) soft nature. 

Lead-Lead collisions at y/s = 2.76 TeV at the LHC provide further information 
about the ridge in heavy ion collisions. With increase in collision energy, the pt 
range of particles that exhibit ridge-type correlations broadens significantly. CMS 
nas observ ed the ridge up to pt — 20 GeV™SI ALICE has also studied the soft 
ridge^EO] anc j their results on centrality dependence are shown on Fig. [3J 

The s tandard explanation of the ridge in HI collisions refers to collective 
flow™. Particularly, angular collimation in cj> can be attributed to the radial 
flow that boosts particles along the radial direction a l ^ . 



54-100% 55-£5% 46-55% 0-5% 




Fig. 1. Soft ridge from STAR The centrality dependence. 



Correlated yield on near-side (IA4>I < 1 ): 
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Fig. 2. Hard ridge from PHOBOS 
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Fig. 3. Soft ridge from ALICE 



CMS 2010, Vs=7TeV 




Fig. 4. Di-hadron correlations from CMS Left: minimal bias data; Right: high multiplicity data. 
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bias event. The right plot shows data on specially selected high multiplicity events. 
The pronounced "ridge" -like structure in the near side (A<fi ~ 0) region extending 
well beyond the jet fragmentation rapidity interval all the way to Ar] = 4 is a 
distinct correlated feature. It has never been observed in pp before, nor can it be 
explained by any of the existing Monte Carlo generators used to simulate hadron 
cascades. 

As was mentioned in the Introduction, only one out of approximately 10 6 events 
leads to extraordinary high multiplicities produced in the final state. It is this set 
of events that exhibits the ridge like correlations. Importantly, the correlations are 
observed only between semi-hard particles with transverse momenta limited to the 
interval lGeV < pt < 3GeV. High density QCD naturally singles out this interval 
as being related to gluon saturation scale Q s ~ 2GeV, typical for such events. 

Following the initial measurement^, further CMS studies have revealed that 
the ridge gets stronger with increase of final multiplicities: the ridge emerges at 
multiplicities of order N ~ 50 and becomes more pronouned at 7-8 times the average 
multiplicity. 

Although multiplicities in the relevant events are large, they are still very far off 
the typical multiplicities in Au — Au collisions at RHIC at 200 GeV. It is unlikely 
that radial flow plays a significant role in these events, and we will take in this 
review the point of view that final state interactions do not change qualitatively 
correlations produced during the collision. If this is the case, the forward angular 
correlations observed by CMS require an independent explanation. 

2.3. Ridge in pPb ©CMS 

Preliminary CMS data on the proton-lead collisions at 5.02 TeV from the LHC test 
run have been released recently^. Analysis similar to one performed by the CMS 
in pp case again revealed the ridge-like structure in the semi-soft momenta interval 
1 < p t < 3 GeV. Special focus given to the ridge's strength (height) as a function of 
the total multiplicity of produced charged hadrons, reveals an approximately linear 
relation (Fig. [5]) , while the overall strength of the ridge is found to be mach more 
pronounced than in the pp case. 

3. Correlations from Theory 
3.1. General considerations 

To explain the origin of the experimentally observed di-hadron correlations one has 
to address two distinct questions: (a) why do the correlations extend over such a 
long rapidity range; and (b) what is the mechanism responsible for correlating the 
observed hadrons within a relatively small angle. 

As has been argued in Ref. ((6]), causality strongly constrains the time at which 
the correlations originate. In order for two hadrons separated by a large rapidity 
interval to be correlated at the time of detection, they must be correlated at very 
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early times: at or just after the collision. As can be seen from the Fig. ([S]), in case 
of HI collisions, the maximal correlation proper time is 

Tlatest correlation Tfreezout 

exp(-Ar//2) (1) 

It follows immediately that the larger the rapidity interval between the two par- 
ticles, the earlier the correlation must be formed. Furthermore, the flatness of the 
ridge in Ar; indicates that all the correlations are produced independently of the 
rapidities of the particles. If this was not the case, the ridge signal would be weak- 
ening for larger rapidity separation between the produced particles, the effect that 
so far has not been observed experimentally. 

In the case of pp collision, the freezout time should be replaced by a (typically 
very short) time scale of parton hadronization. These arguments lead to the con- 
clusion that the long range correlations are present just after the interaction takes 
place, while their production mechanism is practically rapidity independent. 

In QCD there exists a very simple and natural mechanism which can produce 
such long range rapidity correlations - the approximate boost invariance of the wave 




Fig. 5. Multiplicity dependence of the ridge in pPb. 




Fig. 6. Causality structure of correlations 
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function of a very energetic hadron. As we know from variety of experiments as well 
as from the BFKL theory, hadronic wave function at high energy is gluon dominated. 
Since gluons are vector particles, their Lorentz transformation properties are such 
that in a boost invariant state they are homogeneously distributed in rapidity. Thus, 
the observed rapidity independence of the correlations could be merely a reflection 
of the initial conditions of high energy scattering processes as encoded in QCD. 

The angular correlations between produced particles are also to be expected on 
very general grounds at high energy. We will now explain the qualitative picture un- 
derlying these correlations. This is very general and does not rely on specific models 
of high energy evolution and/or hadronic wave function. In the next subsection we 
will provide more calculational details substantiating this qualitative discussion. 

Consider high energy scattering of a hadronic projectile on a stationary target 
in the lab frame. Two projectile gluons with rapidities Y\ and scatter on the 
target and are observed in the final state. As mentioned above, at high enough 
energy the wave function of the projectile is approximately boost invariant. Thus 
the projectile gluon distribution at rapidities Y± and Yi is the same. The two gluons 
scatter exactly on the same target, and thus whatever happens at Y\ also happens 
at Y<z • If for a particular target field configuration a gluon is likely to be produced 
at Y\ at a certain impact parameter, a gluon is also likely to be produced at Y2 at 
the same impact parameter. 

This is especially true in the context of the projectile wave function dominated 
by the large "classical" Weizsacker- Williams field, since in this case fluctuations in 
the wave function are small and the gluon density configuration by configuration is 
almost the same at all rapidities. The hadronic wave function at high energy can 
be computed in QCD using light-cone Hamiltonian perturbation theory™! 

|P) = n|^> = cxp{z J d 2 xb?(x) J dr ] (al a (x,r ] )+at(x,r 1 )) }b(o, a*)\i/>) . (2) 

Here ip is the wave function of valence charges, that determines the distribution 
of fast partons responsible for the color charge density p; B is a Bogolyubov -type 
operator of the soft gluon fields a, and the Weizsacker- Williams (WW) field b is 
given in terms of p via classical Yang-Mills equations of motion. For small (dilute) 
projectile 

&?(*) = - / M*-y)p a (y), /i(*-y) = 7^% s = i.(3) 

I" Jy {x-yr 
For large projectile the WW field is parametrically large b ~ i, while the Bo- 
golyubov operator B produces the fluctuations of the gluon field of order unity. 
Thus for fixed p(x) the gluon density fluctuates very weakly around large average 
value determined by the classical field 

n = (at a) oc b 2 ~ 0{— )■ (n 2 ) - (n) 2 ~ 1 (4) 

Us 

The smallness of the fluctuations is clearly helpful. Although the wave function at 
different rapidities in a boost invariant projectile must be the same, the magnitude 
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of the color field (and therefore the number of gluons) may differ at different values of 
Y for the same configuration of the valence color charge density p, if the fluctuations 
in this wave function are significant. Thus in the same scattering event there may 
be significant differences between particle production at different rapidities. Still, 
although the quasiclassical nature of the state eq.© ensures long range rapidity 
correlations at large values of p, it is not absolutely necessary. Even in the presence 
of considerable fluctuations in the soft gluon wave function, one nevertheless would 
expect positive correlations in rapidity. The only really necessary condition is that 
the density of incoming partons is large enough, so that there is large probability 
to produce more than one particle at a given impact parameter (we will quantify 
what we mean by "given impact parameter" shortly). 

Thus the long range rapidity correlations come practically for free whenever 
the energy is high enough so that the wave function of the incoming hadron is 
approximately boost invariant, and there is very little in the actual dynamics of the 
collision that can affect this feature. 

By almost exactly the same logic we can argue that near-side angular correla- 
tions must be present between the final state gluons. Indeed, if the two gluons hit 
the target at the same impact parameter, their scattering amplitude is determined 
by the same configuration of the target field. Thus if the first gluon is likely to be 
scattered with momentum q, the same is true for the second gluon. One therefore 
expects clear forward correlations for gluons that scatter at the same impact pa- 
rameter. Of course, the two gluons will not scatter always with exactly the same 
momentum transfer even if they hit at exactly the same impact parameter, since 
even a fixed configuration of target fields corresponds to a nontrivial probability 
distribution of momentum transfer. Nevertheless, given that this distribution has a 
maximum at some particular momentum transfer, the angular correlations must be 
very generic. 

The above argument does not invoke the physics of saturation. Nevertheless, 
saturation, if present in the wave function docs affect strongly the gluon produc- 
tion. To better understand the expected properties of the angular correlations in 
the framework of high energy QCD, let us briefly recap our understanding of the 
transverse structure of the hadron in the saturation regime. It is convenient to 
think of the distribution of the (color) electric field configurations in the target. 
The target wave function is characterized by the saturation momentum Q s . The 
saturation momentum plays a dual role in the hadronic wave function. On the one 
hand, it measures the typical magnitude of electric field in the wave function, on 
which gluons scatter. On the other hand it is known that the fiel d comp onents with 
transverse momenta pr < Q s are suppressed in the wave functioj™. This means 
that the electric fields in the target are correlated on the length scale A ~ Qj 1 . 
Thus the saturation momentum doubles up as the inverse of the correlation length 
of target color fields. Typical field configurations in the target can thus be thought 
of having, in the impact parameter plane, a domain like structure of Fig. ([7]). Now 
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consider a projectile parton impinging on one of the domains or "hot spot" of the 
target, with electric field E. While traversing the target field, the parton of charge 
q acquires transverse momentum 



A parton at a different rapidity but with the same charge will pick up exactly the 
same transverse momentum if it scatters on the same " domain" . This of course 
results in positive angular correlation of produced gluons. 

We note that this simple picture also suggests that angular correlations at angle 
(j> and 4> + 7r have equal strength if the hadronic wave function is dominated by 
gluons. Gluons belong to real representation of the gauge group, thus it is equally 
probable to find an incoming gluon with charge q and charge — q in the projectile 
wave function at any rapidity. Suppose for example that on a given configuration 
the color field in the target is in the third direction in the color space Ef = Ei5 a3 , 
while in the incoming projectile the gluon corresponds to the vector potential in the 
second direction A\. One can always write A 2 — — i/2(A — A*), where A = A\ +iA2 
is positively charge with respect to color charge in the third direction, and A~ 
is negatively charged. Thus necessarily equal number of gluons in the incoming 
projectile have opposite sign charges and are kicked in opposite directions while 
scattering on the target. This produces equal strength correlations at angles zero and 
7r. This feature of equal forward and backward correlations was noted in Refs.(|14[ 
I15|) . For quarks which carry fundamental charges, this degeneracy should be absent 
and the presence of quarks in the projectile wave function should skew the angular 
correlation towards zero angle and away from A(f> = tt. 

3.2. A little more hands on. 

The above general discussion provides an intuitive picture for the existence of corre- 
lations. Our next goal is to see how this is reflected in the calculational framework 





(5) 




E 



Fig. 7. Domain structure of chromo-electric fields in the target. 
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based on hdQCD. As explained above, we will base our analysis on expressions 
derived rigorously at high energy for a scattering of a dilute projectile on a dense 
target. We do this for two reasons. One reason is the lack of explicit analytical 
expressions for the case of dense-dense scattering, althou gh a significant adva nce 
in understanding of symmetric collisions has been achievec p9 | 50 | 51 | 52 | 53 | 54 | 55 | anc j 
an elegant numerical approach is being pursued in a series of papers Refs.(|50l |5"TI 
52")) . Also, as we will argue later, the dense-dense scattering approach is not en- 
tirely adequate for study of correlations either, as it does not take into account the 
Pomeron loop effects which in our view is crucial for quantitative understanding of 
correlations. The dense-dilute approach is sufficient for the understanding of basic 
ingredients involved in the calculation of correlations. 

The quantity directly related to the measured correlations, is the double inclusive 
gluon production amplitude. The correlations are computed as variance with respect 
to the single inclusive measurement: 



R 



d 2 N dN dN 



, dN dN 
d 2 kd£ d 2 pdr\ 



d 2 pdr\ d?k d£ d 2 k dt; d 2 p dr\ 

Here p, k are the transverse momenta of the observed gluons, and r), £ are their 
respective rapidities. The experimentally measured correlations are normalized per 
number of produced charged hadrons, R ex P ~ RdN/dt;, however the quantity in 
eq.© has a more straightforward interpretation. 

The expressions for single- and double-inclusive gluon production probabilities 
are known in the literature. We outline the derivation in Appendix A, while here 
we quote the results. The single-inclusive amplitude for production of a gluon with 
transverse momentum k and rapidity r\ reads 

dN 



d 2 kdrj 



nplitudc lor proo 

^561571581591601 



1 



x {p;(x 1 )[SH Xl )~SHz)}[S(x 1 )~S(z)r b p b p (x 1 )} (7) 

whre, as before f(r) = r/r 2 is the Weizsacker- Williams field of a point charge. 

Eq. ([7]) is easy to interpret. The projectile appears in eq.(J7|) as a distribution 
of the color charge density Pp(x). The charge density is a function of transverse 
coordinates x only, since in the high energy limit all the charges are localized at 
x~ = and there is no nontrivial dependence on the longitudinal coordinate. The 
target is described as a distribution of the color fields which enter eq.flTj through 
the eikonal scattering matrix of a single gluon S ab (x) when scattering off the target. 

For a fixed configuration of the projectile charges p a {x\) in the incoming wave 
function, a WW gluon is emitted at a transverse position z with the probability 
amplitude f(x\ — z). This gluon and the parent charge both scatter eikonally on 
a fixed configuration of the target field. The production amplitude of the gluon is 
determined by the eikonal amplitude S(x\) — S(z) which describes the decoherence 
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of the WW gluon from its emitter. The square of this amplitude gives the proba- 
bility of the gluon production in the given event (for a fixed projectile and target 
configuration) . 

This probability then has to be averaged over the projectile and target wave 
functions - the projectile ensemble determining the distribution of p(x) and the 
(independent) target ensemble which determines the distribution of S(x). In the 
framework of CGC, the averaging over the projectile and target ensembles is per- 
formed with the weight functionals W p and W T providing the probability densities 
to find a given projectile and target configurations respectively: 

(a)p, T = J DpW p [p p ] J DSW T [S]a[p p ,S] (8) 

The weight functionals W p and W T are subjects to high energy evolution which 
determines their dependence on rapidity: 

^_ = H RFT W _ (Q) 

Here H RFT is the effective Hamiltonian of QCD at high energies. H RFT has to 
be derived directly from QCD and is the QCD incarnation of the Reggeon field 
theory of the old Reggeon Calculus. The derivation of H RFT at the moment is 
incomplete. For processes involving the scattering of a dilute object on a dense one 
the limits of H RFT are known. The evolution of the wave function of a dense target 
which can scatter on a dilute projectile the relevant approximation is the JIMWLK 
evolution. The complementary form of the evolution applicable for the evolution 
of a dilute target, which participates in the scattering on a dense projectile, is 
the so called KLWMIJ evolution derived in Ref.(|45|). Both of these limits include 
a tripple-Pomeron vertex and effectively resum the fan diagrams of the Pomeron 
Calculus. The JIMWLK and KLWMIJ evolutions are applicable to the same type of 
processes, the asymmetric situation when one of the objects is dense and the other 
one is dlute. The two equations describe the process in different frames: JIMWLK 
ascribes the evolution in energy to boosting of the dense object, while KLWMIJ to 
boosting of the dilute one. A more general form of H RFT which includes Pomeron 
loops is available in Ref. (153[ l54| . however it is not general enough to be applicable 
to dense-dense scattering. 

Eq. ([7]) is not symmetric under the interchange of the projectile and the target, 
which reflects our treatment of the projectile as a dilute object. For a dense target, 
the eikonal factors S are nonlinear functions of the target color charge density, which 
is the result of the resummation of multiple scattering on a dense target field. The 
projectile contribution to eq.Q on the other hand, is quadratic in p p , reflecting the 
fact that the gluon emission amplitude in the projectile wave function is accounted 
for in the linear approximation valid only for small p p . 

The double-inclusive gluon production in the CGC framework has been analyzed 
in Ref.(l6l! 62, 63, 64|. The general expression is rather long and cumbersome, but 
only part of it is relevant to the mechanism we discuss in this review. This is 
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the part of the cross section where the two gluons are produced independently 
from two different showers. Formally this term is leading in the limit of large color 
charge density of the projectile, p p ~ l/g. The relevant part of the double-inclusive 
production cross section is 



Eq. (fTU|) does not include rapidity evolution between the rapidities of the pro- 
duced gluons. It is applicable when the gluons are produced sufficiently close to 
each other in rapidity so that there are not many gluons produced in the rapidity 
interval between them. This is not a real limitation, since the rapidity evolution 
implies that this is the case unless the rapidity difference between the two gluons is 
greater than 0(l/a s ). For a greater rapidity separation one wo uld have to include 
the evolution. The formalism for including the evolution existi^ESl an( j an approx- 
imate numerical study in connection with ridge in HI has been performed^, but 
we will not discuss this more complicated kinematics any further in this review. 

It is straightforward to see that eq. tfTO]) indeed produces angular correlations. 
For a fixed configuration of the projectile sources p p (x) and target fields S(x), the 
function a(k) as a function of momentum has a maximum at some value k = Q. 
Therefore clearly the product in eq. (|10|) is maximal for k = p = Q. The value of 
the vector Q of course differs from one configuration to another, but the fact that 
momenta k and p are parallel to each other does not. Therefore after averaging 
over the ensemble, (<r(fc) &{p)) has maximum at relative zero angle between the two 
momenta. The strength of the maximum of course depends on the detailed nature 
of the field configurations constituting the two ensembles (the projectile and the 
target). 

We stress that even though averaged over all configuration (<j(fc))p ; T must be 
isotropic, there is absolutely no reason for it to be isotropic for any given configu- 
ration. Thus, angular correlations emerge as a result of target / pro jectile averaging 
procedure over isotropic ensembles of anisotropic configurations^. 

The simple formula eq. ([TO)) indeed encodes all the main features of the intuitive 
discussion of the previous section. The correlations originate in the initial state. 
They are clearly rapidity independent within the approximation in which one does 
not include rapidity evolution between the two observed gluons. As noted above, 
this is expected to be valid for rapidity difference r/ — £ < l/a s . 

The correlations are peaked at the forward angle, but, consistent with our pre- 
vious discussion, they also peak at the backward angle with the same strength. 
Clearly, a(k) is a symmetric function of its argument. It is expressed in terms of 
the single gluon production amplitude 



dN 



(er(fc) a(p))p^T + terms subleading in p p 



(10) 



d 2 pd 2 kdrjdt; 



a(k) 



a(k)d* (k) 



(11) 



with 




(12) 
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Since the amplitude a(x) is real in coordinate space (eikonal scattering matrices in 
the eikonal approximation are by definition real), it satisfies a*(k) = a(—k), and 
therefore configuration by configuration 

&(k) = a{-k). (13) 

Thus the leading contribution to the two particle inclusive production probability 
is symmetric under 



cPpcPkdrjdt; ' d 2 pd 2 kdr)dt; 
and must have two degenerate maxima, at relative angles A(f> = 0,7r. 

Eq. fTUl) also suggests that if one includes in the calculation hadrons produced 
from quarks, the degeneracy will be lifted. The quark production amplitude involves 
the eikonal matrix of a produced quark, which, being a fundamental representation 
matrix, is complex. Thus the k — > —k symmetry is broken by the presence of quarks. 

The overall magnitude of the correlation can be estimated from the following 
simple consideration. In order for the two produced gluons to be correlated in the 
final state, they have to be close in the initial state and also scatter off the same 
target field. In other words they must be within the smaller of the two correlation 
length that characterize the projectile and the target wave functions. The target 
and the projectile correlation lengths in the saturation regime are determined by 

i PIT) 

the saturation momenta l p / T \ oc Q s . Thus for correlated production the two 
gluons need to be within the radius \jQ™ ax of each other, where Q™ ax is the larger 
of the two saturation momenta Qf and Qj. On the other hand the total number 
of produced gluons is proportional to the total transverse area o f the smaller of the 
two objects participating in collision. Thus parametricalhJ^^ 

^ ~ (r)max\2q . ' (^) 
s ) u rntn 

Our general discussion provides an intuitive explanation to the correlation phe- 
nomena and also makes it clear that the presence of the correlations does not depend 
on the specifics of the approximation used to estimate ji k ^i p ■ The magnitude of the 
effect however may depend on the approximation quite strongly. In particular what 
one assumes about the projectile and target averaging procedure is very important. 
Numerical calculations of correlations in the CGC approach have used the Gaus- 
sian weight for averaging the projectile and the target charge densities. Although 
the Gaussian weight makes the calculations manageable and explicit, it has certain 
non-generic features which most likely skew the results of the calculation. We now 
address this quite specific, but important point. 

4. Correlations with Gaussian averaging - the large N c 
conundrum. 



For the calculation of the correlated production one needs two types of averages. 
The simpler ones appear already in the single-inclusive gluon production. On the 
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projectile side one needs to know 

( P ;(x)p b p (y) P = J J—S ab {p P (x)p P (y)) P = -2—5 ab <S> P (x,y) (16) 

$ P is related to the unintegrated gluon distribution of the projectile via Fourier 
transform: 



$ p {k 1 ,k 2 ) = [ e- lfclMfc ^$ p (x ;2/ ) 

Jx,y 



(17) 



The color structure of eq. (fT6"| follows directly from the color singlet nature of the 
projectile wave function. Assuming translational invariance in the transverse plane, 
one has 



<S> P {x,y) = <$> p (x - y); $ p (fci, k 2 ) = p (fci) 5 2 {k 1 + k 2 ) 



(18) 



This approximation is commonly used in the literature. It should be a good approx- 
imation whenever the hadron is large enough to contain several domains of the size 
Q~ s \ 

On the target side the relevant average is 



{[S{x)sHy)] ab ) T = 



l 



-8 ab {tr[S{x)S\y)] 



(19) 



Here tr[S(x)S^(y)] is the forward scattering matrix of an adjoint (gluonic) dipole. 
Expressed in terms of its Fourier transform 

1 



T(k u k 2 ) = i 
the single inclusive probability reads 



-ik\x—ikiy 



(tr[S(x)S\y)]h 



dN 
cPkd-q 

with the vertex 



— / K 2 ^2{k 7 p,q)$ P (p,q)TY-n{-k-p,k-q) 
n Jp,q 



K 2 ^ 2 (k,p, q) 



qi k{ 




ki 


_ Pi_ 


q 2 k 2 




k 2 


p 2 



(20) 



(21) 



(22) 



Here Y is the total rapidity interval between the target and projectile (total energy 
of the process) while 77 is the rapidity of the observed gluon. If either the target 
or the projectile are translationally invariant, Eq. (|2ip reduces to the standard k t 
factorized formula for DIS I5S1. 

The other type of averages involves the fourth power of the gluon eikonal am- 
plitude in the target: 

([SHx)S(z)} ab [SHy)S(u)} cd ) T (23) 
and similarly fourth power of the charge density in the projectile 

(p a p (x)p h p (x) P ;(y)p d p (y)) P . (24) 
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In the Gaussian approximation^!^ one assumes that W p (and similarly for 
W T , see below) is a Gaussian functional in p at any rapidity. This leads to the 
factorized correlation functions: 

(p a p (x)p b p (x) P ;(y)p d p (y))p = [(p a P (x)p b P (*))p(p c P (y)p d P (y))r + 

+ (pl{x)p c P {y))p(p P {x) P d P (y))p + (p»p d P mp(p P {y)p b P {x))p] 

= m ~ 1)2 [8 a *S cd * P (x,x)* P (y,y)+ 6 ac 6 bd <i> p (x,y)<l> p (x,y)+ 
+ 6 ad 5 bc $ p (x,y)$ p (y,x)] (25) 

Once we assume eq.(|25|) for the projectile averages, only three target structures 
enter the calculation: 

x (tr[(SHx) S\z))(S(x) - S(z))]tr[(S\y) S*(w))(S{y) - S(w))}) T + 
+$ p (x - y)$ p (x - y) x 

x (tr[(S\x) - Si(z))(S(x) - S(z))(S\y) S^w))(S(y) - S(w))]) T + 
+$ p (x - y)$ p (x - y) x 

x<tr[(St(„) - S\w))(S(y) - S{w)){S^{x) - S\z))(S{x) S{z))]) T 

(26) 

This expression contains two types of target matrix elements: the two-dipole 
matrix element (tr[S^ S]tr[S^ S]) T and a quadrupole (tr[S* SS^ S]) T . 

In principle one is interested in dependence of the double gluon production on 
both, the rapidity of the observed gluons and on the total energy of the process. To 
study these one has to evolve the projectile averages to rapidity -q and the target 
averages to rapidity Y — r\ with H RFT . Current numerical implementations adopt 
the simplified approach, whereby the projectile averages are evolved by the BK 
equation, while for the target averages one adopts an ansatze inspired by the limit 
of dilute target. In the dilute target limit one can expand the eikonal factor S to 
leading order in the target color field 

S(x) = l + iT a a a T {x) (27) 

where af. is the color fields created by target charges, V 2 ar = pr- One further 
assumes Gaussian weight for the target color charge densities, or equivalently for 
the color fields oct 

{a a T {x)a b T {x)a c T {y)a d T {y)) T = (jy2 \ X (28) 
6 ab 5 cd $ T (x, x) $ T (y, y) + 5 ac S bd $ T (x, y) $ T (x, y) + 6 ad 6 bc $ T (x, y) $ T (y, x) 



where 

<!> ' /• x\ = 

V 2V ' V 



<S> T (x,x) = — ( x -v)^(x- v)<S> T (u,v) (29) 
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The rapidity dependence of the averages is then given by substituting for $>Y-n the 
solution of BK equation. 

This procedure amounts to the following ansatz for the target averages at all 
rapidities: 

({S{x)S^y)Y\S{u)S\v)y d )=5 a H cd — ±^{tr[S{x)S\y)]){tr[S{u)S^v)]) 

+ {T l TT\TT>y d N2{ ^ 2 i)2 {tr[S{x)S\u) - l]){tr[S{y)S\v) - 1]) 

+ (T i Tiy b (TiT i y d N * _ (tr[S(x)&(v) - l]){tr[S(u)S\y) - 1]) (30) 

Substituting eq. (1301) into eq. (|26p gives an explicit expression for the double inclusive 
probability in terms of solutions of the BK equation. The structure of the final result 
is summarized in a self-explanatory Fig. (j4]), borrowed from Ref. (|70|). 

We will not write out all the terms explicitly, but rather stress one very important 
point. Namely, the correlated part of production probability in this approach is 
subleading in N c . This is easy to see, since the leading N c contribution to the 
doubly inclusive probability in the Gaussian approximation is given by 

d 2 vd 2 kdjld£ ~/ >< (31) 

ll jjuj n. u,//u,^ Qaussian J x,x,y,y.z,z,w,w 

x e^-^$ p (x - x)(tr[{S^(x) - S^(z)){S{x) - S(z))]) T x 

x M w -^$ p {y - y)(tr[{S\y) - S\w)){S{y) - S{w))]) T 
dN dN 
d?pdrj d 2 kd£ 

This contribution therefore cancells in the experimental ratio R, and only the terms 
suppressed by l/N 2 contribute. This leading N c cancellation is specific to the Gaus- 
sian approximation, and is not generic. We will discuss it in more detail in the rest 
of this review. 




Fig. 8. Correlated two-gluon production within the Gaussian approximation 
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Gaussian approximation is attractive due to its relative simplicity: it allows to 
compute particle correlations using the known blocks computed earlier and fitted to 
other data. The first application of this approa ch gave reasonable agreement with 
the CMS pp datJSl. More recent paperJZD^l include the " 'nonclassical" ' terms 
that we have not explicitly written in eq. (|10p to successfully fit pp data and give 
a reasonable prediction for pPb double production not just in the ridge kinematics 
but at all angles. 

The theoretical status of this approach in the present context is however rather 
questionable. For one thing, as mentioned above we expect non-vanishing leading 
N c contributions to the correlated production. The other aspect is, that even in the 
subleading order in N c , the Gaussian approximation does not give a generic result. 
In particular as was pointed out in Ref. (|70| . the structure of the averages eq. (|25|) is 
not compatible with the BFKL and BK evolution. It was suggested in Ref. [70] that 
additional terms of the form 

Hp> h yPuPi) = P^[8A]M(x, y, u, v) + P t a c b d [8A]M(x, u, y, v) + Pf cd [8A]M(x, v, u, y) 

(32) 

with the projector P[8A] denned in eq.([34|). should be included in the averag- 
ing procedure. One can in fact perform the complete analysis of the in dep endent 
"condensates" that can appear in the projectile (target) wave function-^. Given 
a color singlet projectile (target), the most general expression for the matrix ele- 
ment is expandable as a sum over a full set of irreducible representations to which 
8 <g) 8 = 1 + 8a + 8s + 10 + 10 + 27 + R 7 can be decomposed (see e.g. Ref. (|74|): 

(p;(x)p h p (x)p c p (y)p d p (y))p = £ j-i#M A( r] (x,x,y,y) 

r 

Af r] (x,x,y,y) ee (P^> n [r] p l {x)p> ' {x)p m {y)p n {y))p (33) 

Here P[r] is a projector on representation r of dimension d[r] which can arise from 
multiplying two adjoint representations. For example 

P[l] = J^-jS ab S cd ; d[l] = 1; P[8 A ] = ±-f aH f icd ; 

P[8 S ] = -Ld aM d zcd ; d[8] = iV c 2 - 1 (34) 

^ * c 

The explicit form of the other projectors can be found in Ref. (74). In the dilute 
projectil e limit the rapidity evolution of these condensates is governed by the BKP 
equatiorJ^EH, a generalized version of the BFKL equation satisfied by reggeons 
containing more than two reggeized gluons. 

With a similar decomposition for the target matrix elements 

([S^(x)S(z)] ab [S^(y)S(u)} cd ) T = £ ~^ P cd\- r } ^ir]( x > z >yi u ) 

r 

AUx,z,y,u) ^ (P% n [r] [5t(o:)5(z)]^[5t(y)5(u)] m ") T (35) 
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we write 
dN 



cPpd 2 kdr)d£ ^— ' d r 

[Af r] (x, x, y, y) - Af r] (x, x,y,w) - Af r] (x, x,w,y) + Af r] (x, x,w,w) 

- Ajyj (x, z, y, y) + Af r] (x, z, y, w) + Af r] (x, z, w, y) - Af r] (x,z,w,w) 

- Ajyj (z, x, y, y) + Af r] (z, x, y, w) + Af r] (z, x, w, y) - Af r] (z, x, w, w) 
+ a h ( z , z, y, y) - Af r] (z, z, y, w) - Af r] (z, z, w, y) + Af r] (z, z, w, w)} 

(36) 

In Eq. ([56]) we have not indicated the WW field factors, the Fourier phases and 
the integrals over the transverse coordinate for the sake of brevity. Eq. (|36|) does 
not presuppose Gaussian factorization, but if it is assumed, the previous results are 
recovered. The leading N c contribution comes from the singlet projector (r = 1). 
The rest of the terms are N c suppressed and provide the general expression for 
the 1 /N c corrections beyond the Gaussian approximation. To calculate the conden- 
sate one would need to have a model at initial rapidity and subsequently evolve it 
with rapidity using H RFT , or, if suitable, its JIMWLK limit. A model for initial 
conditions that goes beyond simple Gaussian was discussed in Refs. ([771 175]) . 

The rapidity dependence of the condensates is a separate interesting question. 
The evolution of the quadrupole operator (which is a linear combina tion of t he 
condensates defined above^2l has been a subject of active recent studie J 79 l 8 ° ISli. It 
is governed by the JIMWLK equation. Given the fact that the quadrupole in R 
emerges only at order 1/N%, its rapidity evolution can be studied at the leading 
order in N c . Such an evolution was first derived in Ref. (61) and then later in Refs. 



It is also important to remember, that in order to obtain the complete 0(1/ N^) 
result one not only should include direct contributions of the quadrupole and other 
condensates to the observable, but also take into a ccount the 0(1/ N 2 ) effects in the 
evolution of the leading two-dipole operatoJ^EH The latter involves quadrupoles, 
sextupoles" 1 ^ and, likely also higher order multipoles given the fact that it is the 
adjoint dipoles that appear in the expression for R. 

Generally speaking, beyond the Gaussian approximation, none of the conden- 
sates A r ; factorize. In particular there is no a priori reason to assume that the 



factorization holds for the leading condensate 

Af x] (x,x;y,y) ^ j^- ^ (x - x) (y - y) 



Thus one does expect nonfactorizable contributions, and therefore correlated hadron 
production in the large N c limit. 

One might think that Gaussian factorization in the large N c limit is natural 
due to presence of large number of degrees of freedom. However, even though the 
number of degrees of freedom in large N c limit is large, the theory has legitimate 
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states which contain small number of particles. A color dipole is an example of 
such state. It is a superposition of many states (different color orientations) of two 
particles, rather than a state with many particles. This type of state maintains high 
level of color coherence in the large N c limit, and is not subject to the central limit 
theorem. Fluctuations in density in a state like this can be large even at large N c , 
and large fluctuations in the ensemble break factorization of correlation functions. 

A similar question was considered a while ago in Refs.([ 55 1 l57 | in connection with 
factorization of dipole densities in the dipole modeP^^Mll As shown in Refs. 11861 
I87p within the dipole model (which is defined entirely within the large N c limftP^) 
the product of two densities does not factorize. When evolved from a single dipole 
initial state, the connected correlator of dipole densities of small dipoles (of size r, 
much smaller than that of the initial dipole ro) as a function of the distance b r 
between the two, behave as 

(n(r, 0)n(r, b)) - (n(r, 0)) (n(r, b)) = (n(r, 0)) (n(r, 6)> (37) 

where A > is a number. This is not quite the regime of interest to us, since we 
need to consider the situation when both dipoles are at the same impact parameter. 
Refs. (86, 87) also did not consider angular correlations. However it demonstrates as 
a matter of principle, that large N c limit does not eliminate correlations. 

Gaussian averaging procedure described above is much more restrictive than 
large N c limit, and it is very likely that it tends to underestimate correlations. 

5. The dipole model at large N c . 

The discussion in the rest of this review will be devoted to attempts to understand 
the correlations which are leading in the large N c limit. The leading N c contribution 
is proportional to 

Aft (tr[S( Xl ) &(x 2 )] tr[S(x 3 ) S^(x 4 )]) T 

As discussed in detail above, the region of the phase space responsible for corre- 
lations is when all four points Xi lie within the correlation length 1/Q S of each other. 
In this regime one does not expect Gaussian like factorization of averages. Thus one 
needs a model of the target that encodes such correlations. It is also important to 
understand how the correlations evolve with energy. The purpose of this section is 
to describe an initial study of these questions performed in Ref. (|9ip . 

We study the rapidity evolution of an initial ensemble which encodes nontrivial 
correlation, within the projectile dipole model. This model is a leading iV c approx- 
imation, and therefore fits our purpose well. Although the actual observable of 
interest is the correlated gluon production, which is directly related to the average 
of two adjoint dipoles, for the sake of simplicity we will study correlations of fun- 
damental dipoles. The fundamental dipoles s(x, y) are the basic degrees of freedom 
of the dipole model. 
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Our strategy is to choose an ensemble Wq^s] of initial configurations s(x,y), 
which contains nontrivial angular correlations. Each configuration of the ensemble 
is evolved independently according to the BK equatior 

ds(x,y) _ a s f (x - y) 2 



, v ., , , m ^ 1 s ( x > y) - s ( a: ' z ) s (y> z ) ] ( 38 ) 

aY 2ir J z (x — zy (z — y) £ 

with a - the 'tHooft coupling, which is finite at infinite N c . The correlations at the 
final rapidity are the n calculated by averaging the correlator over the ensemble of 
solutions sy(x,yf^. 

DsW$[s]s(x,y)s(u,v) = J DsW£[s]s Y (x,y)s Y (u,v) (39) 

where s Y (x, y) is the solution of the BK equation (|38p with initial condition s(x, y). 
In order to study angular correlations we have to allow the individual members of 
the initial ensemble to be anisotropic. The rotational invariance is then restored 
by averaging over the whole ensemble and not configuration by configuration. In 
Eq. (|38|) s(x,y) is the dipole scattering amplitude for a single configuration of the 
target field, rather than a target average. 

Like in most BK studies, we take the target to be translationally invariant, 
configuration by configuration. That is, we take the dipole scattering amplitude for 
all configurations of the ensemble not to depend on impact parameter, s(x, y) = s(r), 
with r = x — y. The reason for this choice is solely the simplicity of implementation. 
One certainly should study the impact parameter dependent configurations. In the 
context of the present work we view the impact parameter independence as a proxy 
to studying dipoles that scatter at the same impact parameter. Intuitively one 
expects impact parameter variation not to change the behavior of dipoles scattering 
at the same impact parameter too strongly. This question however has to be studied 
separately. 

We start the evolution with a simple one parameter ensemble at some initial 
rapidity Y$: 

s(Y ,r) = exp{-(rE) 2 } = exp{ — r 2 E 2 cos 2 (0 — 5)}; (40) 

The parameter E represents the chromo electric field in a given target configuration. 
The field E is not colored and thus cannot be associated directly with the color field. 
Instead a better way to think about it is as representing color averaged moments 
over the (color invariant) target ensemble, e.g. EiEj ^< EfEj >. The magnitude 
of E plays the very same role as Q s . 

Here 9 and S are the angles of the dipole moment and the field relative to an 
arbitrary axis. The absolute value of the field E was fixed in this study. The angle 
S was taken to be a random variable, so that the weight function for the target 
ensemble is a constant, Wy Q (S) = l/2ir. 

The ensemble we work with, albeit very simple, differs significantly from previous 
numerous studies of the BK equation. In these studies the initial conditions are 
parameterized as a function of the dipole size only, without the angular dependence. 
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This includes papers which do allow for fluctuations in the initial ensemble. For 
example Refs.(92, 93) allow for the fluctuations of the strength of the field E but 
without angular dependence. Such an ensemble, of course, is not suitable for a study 
of angular correlations. Eq. (|40l) on the other hand provides a two-dimensional initial 
data set. 

Fig. (J9j) displays the initial condition plotted for the imaginary part of the dipole 
scattering amplitude N(f,S = 0) = 1 — s(r). It is quite obvious that the initial 
conditions (j40|) with the uniform distribution of 8 generate angular correlations 
between two dipoles. In particular for small dipoles (r, r' <C Qj 1 ) 



Each configuration of the ensemble is evolved independently according to the 
BK equation. The correlations at the final rapidity are then calculated by averaging 
the correlator over the ensemble of initial conditions according to (|39p. 

We observed in Ref. (j9"l"T) that for any given initial condition, the angular depen- 
dence of N(f) is washed away exponentially quickly with rapidity. As a consequence, 
the initial correlations also quickly disappear with rapidity. In the range of rapidities 
studied in Ref.(91) this fall off is approximately exponential 



This result is somewhat surprising in view of our previous arguments. We would 
expect a similar fall off for large dipole sizes, namely those sizes that cross the 
saturation boundary during the evolution r > Qj 1 (Y), but not for small size dipoles 
which stay inside the saturation radius for all rapidities studied. The fall off observed 
in Ref.(91) on the other hand, is rather uniform in dipole size and takes place also 
for very small dipoles. 

Thus it appears that, even if the initial ensemble contains nontrivial leading 
N c correlations, after a short rapidity evolution the correlations disappear. This is 
consistent with recent numerical simulations and analytic ansatze for solutions of 



(s(r)s(r'))s ~ rV 2 +2(r-r*) 2 . 



I? 




A ~ 0.6 



(41) 



N 



0.8t 

o .a 

0.2 

o.c 




Fig. 9. Initial conditions. Left: profile of N as a function of r and 8. Right: the same in polar 
coordinates; the "blackness" of the target is shown with respect to the dipole's orientation. 
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the JIMWLK equation. T he solution seems to converge to ensembles with pretty 
much Gaussians properties2^25129J80| 

At first sight this is inconsistent with the results of Refs.(86, 87), as quoted 
above. However the physics of this decrease, as elucidated in Ref . ([91]) . is fairly 
simple and is not in contradiction with Refs. (|86l |87|) . The calculation of Ref. (|9Tj) is 
performed in the framework of the projectile dipole model, while that of Refs. (|861 
|8"T|) in the target dipole model. Technically the projectile dipole model assumes that 
the projectile is dilute and contains a small number of dipoles, while the target is 
a dense system, while the target dipole model reverses the roles of the two objects. 
In terms of the target evolution, the projectile dipole model corresponds to the 
large N c limit of the JIMWLK evolution, while the target dipole model to that of 
KLWMIJ. 

As discussed in Ref. (46) the JIMWLK and KLWMIJ evolutions are very different 
as far as the target wave function is concerned. KLWMIJ is the normal perturbative 
evolution of the target wave function. It is well known that in such evolution the 
number of gluons grows exponentially fast. The target gluon density, which in this 
regime satisfies the BFKL equation, depends on the total rapidity as 

g(p,Y) cx e ca ° Y (42) 

where c is a number of order one and Y is the total rapidity by which the target 
wave function has been boosted from " rest" . The exponential dependence also holds 
for the differential gluon density in the wave function at any rapidity r\ < Y , namely 
(we suppress the momentum dependence in the prefactor for simplicity) 

ife^e^-,). (43) 
drj 

Gluons in the wave function which are separated in rapidity by no more than 
5rj ~ 0(1/ a s ) are expected to be correlated. Since the gluon density in the wave 
function exponentially grows with rapidity, any projectile that probes such a wave 
function effectively feels only the gluons in the last rapidity "bin". Thus such a 
probe is sensitive to any correlation that exists between the softest gluons in the 
wave function. 

The rapidity dependence of gluon density generated by the JIMWLK evolution 
on the other hand, is very different. As discussed in Ref. (46), while the probability to 
emit an additional gluon in the dilute regime is proportional to the number of gluons 
in the wave function, in the dense regime this probability approaches a constant. 
The evolution with constant probability of emission generates gluon density which 
is uniformly distributed in rapidity. This can be thought of as a random walk in 
color space as in Ref. (|4"5jl . So now one has 

4-9M = C. (44) 
drj 

with C a function of transverse momentum, but not of rapidity. It should still be 
true, like in the KLWMIJ case, that gluons separated by a small rapidity interval are 
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correlated, while the correlation disappears for gluons at very different rapidities. 
However now, as opposed to KLWMIJ case, if one scatters any projectile on such a 
target, the projectile will sample gluons at all rapidities equally. Thus, for example, 
if the projectile consists of two dipoles, the two dipoles will most likely scatter 
on color field components (gluons) with very different rapidities. Since such fields 
are not correlated, the two dipoles will scatter independently and the two dipole 
scattering amplitude will not exhibit any correlations. 

Another way of understanding the difference between the nature of JIMWLK 
and KLWMIJ evolutions is the following. In KLWMIJ evolution one starts initially 
with the target wave function which contains a small number of gluons. The prob- 
ability to emit an additional gluon in one step of the evolution is small (0(a s )), 
however if a gluon is emitted, the end configuration strongly differs from the ini- 
tial one, since the number of gluons have changed by a factor of order unity. Such 
evolution thus generates a very rough ensemble of target field (gluon) configura- 
tions. In this ensemble configurations with very different properties are present, 
albeit with small weight. Such an ensemble exhibits large fluctuations in a variety 
of observables. 

On the other hand, in JIMWLK evolution the probability to emit an extra 
gluon is large - of order unity. However emission of an extra gluon produces a new 
configuration which is hardly distinct from the existing one, since the one extra 
gluon is produced on the background of — gluons which already exist in the wave 
function. Thus the JIMWLK ensemble contains many configurations, but all these 
configurations have very similar properties. The fluctuations in observables in such 
an ensemble are very small. 

Thus it is reasonable to expect that KLWMIJ evolution produces (and pre- 
serves) correlations (and fluctuations) between scattering amplitudes of two projec- 
tile dipoles at leading order in N c . These are indeed the correlations discussed in 
Ref. (|86l |87|) . On the other hand JIMWLK evolution must lead to disappearance of 
any correlations initially present in the target ensemble, and this is exhibited by 
the results of Ref . (HJ) . 

Given that the two approximation to high energy evolution lead to such qual- 
itatively different answers for the observable we are interested in, naturally one 
should ask which one of them, if any, should be used for quantitative investigations. 
Naively one might think, that since we are interested in high multiplicity situation, 
the target is dense and JIMWLK evolution is more suitable. This however is not 
the case. We expect the correlations to arise due to scattering at close values of 
impact parameter. The transverse distances in question should be smaller or of the 
order of the saturation radius Qj l . This means that the scattering occurs on the 
components of the target color field with transverse momenta p > Q s . However 
at these momenta the target wave function by definition is still dilute. Even for a 
dense target JIMWLK evolution is not appropriate for all wavelengths, but only for 
those greater than the saturation radius. At large momenta p > Q s the appropriate 
evolution is KLWMIJ and not JIMWLK. 
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Thus to be able to describe fluctuations one needs to be able to evolve the high 
momentum modes according to KLWMIJ while low momentum modes according 
to JIMWLK. In fact, the situation is even more complicated. One expects that the 
largest contribution to correlations comes from the modes with p of order Q s , since 
in the relative distances of order Qj 1 should dominate the integral over impact pa- 
rameter. For these modes however, the Pomeron merging and splitting (JIMWLK 
and KLWMIJ contributions) are of equal importance. In other words, this momen- 
tum range is sensitive at the leading order to the Pomeron loop contributions. We 
thus conclude that Pomeron loops need to be included for proper treatment of this 
question. 

The impact parameter independence of the distribution considered in Ref. (|91j) 
obviously misses an important part of dynamics, as it largely based on hot spot 
picture of the transverse plane. The full picture should include the electric field 
E = E(b) fluctuating both in magnitude and direction throughout the transverse 
plane (Fig. [7]) . This study is currently in progress^'. 



6. Double Pomeron exchange model 

As mentioned in the introduction, there is another approach to angular correlations 
in the recent literature^, that seems to be rather close to the one discussed here. 
In this section we wish to discuss it somewhat further, and also point out that 
superficial similarities notwithstanding, the physics described by the approach of 
Ref. (|Tr?|) is quite different. We will not follow Ref.([TF|) verbatim, but rather will give 
our interpretation of the approach presented there. 

Let us examine more carefully our original expression eq. (|10[) . For the purpose 
of the present argument we will assume that the production is dominated by the 
part of the phase space where the charge densities p{x{) are far enough from the 
emitted gluons, so that the gluons necessarily decohere and materialize in the final 
state once they scatter. This is essentially partonic picture of the projectile, where 
the gluons in the incoming wave function are assumed to carry only very small 
transverse momentum, and all the momentum in the final state comes from the 
exchange with the target!^. Although this contribution does not account for all the 
gluon production, it is likely to dominate the correlated production at <f> ~ 0. With 
this assumption we can simplify our expressions by dropping the eikonal factors 
associated with the scattering of the valence color charge densities. The expression 
for double inclusive gluon production then becomes 

dN a 2 f d 2 l d 2 l d 2 Q 

(tr(fc) a(p)) PtT = -2 / 7^17^7^ x 



d 2 V d 2 kdj]di x v ; ir 2 J (2tt) 2 (2tt) 2 (2tt) 2 



k 2 {l-Q) 2 J \k 2 (l + Q) 



Pj h + Qj \ ( I'' 1 



P 2 (1 + Q) 2 J\p 2 (l-Q) 



x (n(l, Q)h(l, -Q)) P {f(k - I, -Q)f(p - I, Q)) T (45) 
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with 

n(l,Q) = J d 2 xd 2 ye a{x - y)+tQ(x+y) p a (x)p a (y)- 
f(k,Q) = ( d 2 xd 2 ye lk{x - y)+lQ{x+y) TY[S{x)S\y)] (46) 



Here n (T) are operators on the projectile (target) Hilbert spaces such that 
<S> p (-l-Q,-l + Q) = (n(l,Q)) P ; 
T(-l-Q,-l + Q) = w ±- J (f(l,Q)) T (47) 

In the above equations we have assumed that the projectile and the target 
are both described by translationally invariant ensembles. Although we have not 
indicated this explicitly, both h and T depend on rapidity, so that strictly speaking 
one should write n v , Ty- V - Note, that the physical meaning of Q is not that of 
momentum transferred in the scattering, but rather that of the initial transverse 
momentum of the dipole with fixed size r. Thus T(Q) has the meaning of forward 
scattering amplitude of the dipole with transverse momentum Q. 

For single inclusive gluon emission amplitude we again quote eq. ([7]) 



dN a f d 2 l 



d 2 kdrj 7r J (2-7T 



k{ \ ( ™i li 

k 2 ^! 2 \¥ ~ P 



(n(l,Q = 0)) P (T(k-l,-Q = 0)) T 

(48) 

The double inclusive cross section can be written in terms of a natural generalization 
of eq. (|4"8"]) . With minor abuse of notations we introduce an operator 



dN(Q) _ a f d 2 l 



d 2 kdr\ 7r J (27T 



ki 1% Qi \ ( k{ li ~\~ Q 



k 2 (l-Q) 2 J\k 2 (l + Q) 



n(l,Q)T(k-l,-Q) 
(49) 

keeping in mind that eq. (|49l) is defined configuration by configuration for fixed values 
of projectile density and target field. This quantity (up to kinematical factors) can 
be though of as a Pomeron extending from the projectile to target rapidity which 
is cut at the rapidity where the final state gluon is measured. 

The translational invariance of the target and projectile therefore restricts 
the contributions to the single inclusive cross section to zero momentum transfer 
Pomeron. For double inclusive gluon however, only the total momentum transfer 
between the two Pomerons has to vanish, while each Pomeron individually can 
transfer a nonvanishing momentum Q. Thus a priori, there are two distinct sources 
of nonfactorizability of the double inclusive cross section, and thus of correlations. 
First, is possible nonfactorizability of the Q = amplitudes, and second direct 
contribution of the Q ^ terms. 

In the calculation described in the previous sectiorJ^, by limiting ourselves to 
the target ensemble translationally invariant configuration by configuration we have 
limited ourselves to the Q = contribution, and found that the initial correlations 
fall off exponentially with rapidity. On the other hand Ref. (fT6f considers the effect 
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of Q ^ contributions within the framework of the BFKL evolution and finds 
correlations which survive at arbitrary rapidity. 

6.1. Disappearance of correlations at Q — 0. 

Within the BFKL framework at leading N c initial correlations possibly present in 
the Q = channel, disappear with rapidity. Consider for definiteness the target 
averages. In the large N c limit, the solution of the BFKL (or rather BKP) equation 
for (T(k, Q)T(p, —Q))t can be written in terms of the eigenfunctions of the BFKL 
equation for the single Pomeron. Recall that the complete basis of solutions of the 

[no] 

BFKL equation is 

( (x-y) / (x-y)* \^ +iv 

= {(*-pHv-p)J {( x -py(y-py) (50) 

where the complex coordinate x is defined as x = x\ + 1x2 The corresponding 
eigenvalues do not depend on p 

u n , v = 2a s x(n,v), (51) 

where 

X(n, v) = *(1) - \ * (1±M + i„) - \* (i±M - <„) , (52) 
These wave functions satisfy the completeness relation 

16 ( v 2 + 4 



^ pOO p 

(2ir)H(x-x)5(y-y) = £ / dv / d l p 

m m J — OO J 



71— — OO 



x - y\ 2 \x - y\ 2 



x ^( x - p ,y- p )^«*(x-p,y-p) (53) 
The momentum space representation of the BFKL eigenfunctions 

Vw, P (fc,Q) - / d 2 (x~y) J d 2 {x + y)e lk < x - y)+lQ -^iP n ^ p (x,y). (54) 

For translationally invariant solutions at Q = the basis is 

i> n ,u(k) = J d 2 (x-y)e ik( - x -^ J d 2 p i> n ,v, P {x,y) (55) 

At large N c we have 

(T(k, Q = 0)T(p, Q = 0))t = J dv\dv 2 ?.n®{vi,V2,n)\l) n , Vl {k)il>- n , V2 {p) (56) 

S n $(n)V„(A;)'0-n(p) 

where we have assumed rotational invariance of the target ensemble. We have also 
assumed for simplicity that the integral over v\^. is dominated by the saddle point 
v\ji = 0, since this corresponds to highest eigenvalue at every n. This assumption is 
not strictly valid for arbitrary transverse momenta. In particular one expects that 
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at high k,p the saddle point will depend on the values of momenta^, but we do 
not expect this dependence to change our argument. Due to the presence of n 7^ 
terms this sum contains angular correlations between the two densities. At the 
initial rapidity these correlations, as encoded in our initial ensemble determine the 
coefficient function $(n 7^ 0). In particular, the initial distribution ea. tjUI)) contains 
all even angular momenta n. Within BFKL evolution each term in the sum in eq.(|56p 
evolves with rapidity with its own exponential 

$y(n)oce 2w " (y - yo) $y (n) (57) 

The eigenvalue uq is the standard BFKL intercept and is positive, while all other 
eigenvalues are negative. For v = 

w„=o = 41n2fi s ; w n=2 = 4(ln2 - l)a s (58) 

Thus the re/0 contributions decay exponentially and after short rapidity interval 
the n = component totally dominates the sum 

(f Y (k,Q = 0)f Y (p,Q = 0)) T -^y^ ~ ^(k)^(p) (59) 

For n — the function ip(k) does not depend on the direction of momentum k and 
therefore eq. (|59[) does not contain any angular correlations. The result quoted in 
(|41[) can be indeed identified with Lo n =2 evaluated at a reasonable value of a s ■ This 
is indeed the origin of the exponential disappearance of angular correlations in Ref. 
flS]). 



6.2. The Q ^ contribution. 

Ref. (16) considers the general case of Q 7^ and n = 0. In this case one has 
{T(k,Q)f(p,-Q)) n =a= / (fpx / d 2 p 2 \ $ UltV2 (pi, P2)^u 1 , Pl (k, Q)^u 2 , P2 (p, -Q) 



^ d*pi I d"p 2 $(puP2)i> Pl (k,Q)i> P2 (p,-Q) (60) 
where 

= / d 2 (x-y) [ ^(x + y)e a -^^- Ss ^ 1 l> n=0 , VtP (x,y) (61) 



and in the second line of (|60l) we have again assumed that the v integral is dominated 
by the saddle point = (tp p = Vv=o,p)- In general eq. ([60l) contains angular 
correlations. Since ip depends on the impact parameter b = x+y only in combination 
with p, one can shift b by p in the Fourier integral to obtain 

(f (k, Q)f(p, -Q)) n=Q = J d 2 p 1 d 2 p 2 e^~^^( Pl ,p 2 )^ p=0 (k, Q)^ P=0 (p, -Q) 
= $(Q)V P =o(fc, Q)^=o(p, -Q) (62) 

Note that 

ip P =o(k,Q)=iP(k 2 ,Q 2 ,k-Q) (63) 
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This expression correlates the direction of k with that of p, through the dependence 
on Q. This is true even for the rotationally invariant target ensemble, which in the 
present context means &(Q) — &(Q 2 ). 

In particuar Ref.(|16p assumes for illustrative purposes 

dN(Q) dN(Q = 0) g ~ a dN 
dk 2 dr) dk 2 dr] dk 2 drj 

Then the correlated part of the double gluon production is 

d 2 N f , Q o, ,r ^o,- *>o, dN dN 



dk 2 dr/d 2 pd£ i 



correlated 



dN dN 
x dk 2 dr\ dp 2 dt;' 

Thus one finds correlation, which conforming to the earlier discussion has a peak 
when the two momenta are either parallel or antiparallel. Since the expression in 
eq. (|62|) involves only the BFKL eigenfunctions with n = 0, it has the leading energy 
dependence and in the dilute regime grows with energy with the double Pomeron 
intercept. Thus the angular correlations encoded in eq.(|62j) survive at high energy. 

So it looks like the result of Ref. ([T6|) is very much in line with the qualitative 
picture we outlined in the first sections of this review. However, the situation is not 
quite that simple. Consider a single Pomeron amplitude at fixed impact parameter 



T(r,6) = / d 2 P P{p) = / d 2 P P(p- b) (66) 

J \r-b-p\\f+b + p\ J \r-p\\r + p\ 

Now suppose that the function P has a maximum for some value of its argument, 
and for the sake of definiteness let us choose this value zero. Then the integral is 
dominated by p = b, and we get 

a - r 2 

T(f, b) oc = (67) 

\r- b\\r + b\ 

For small values of the dipole size (smaller than the impact parameter b) 

2 

f(r,6)cx^ (68) 

This is just the rotationally invariant scattering amplitude of a small dipole on a 
density profile that decreases as 1/b 2 . Intuitively a dipole will scatter this way on 
an ensemble of fields, whose magnitude decreases according to Coulomb law away 
from a central charge at b — 0. The ensemble itself is rotationally invariant, which 
leads to an isotropic scattering amplitude. On the other hand for r 2 ^> b 2 we have 

T(r, b) — const (69) 

So a dipole scattering in the center of the distribution again has an isotropic scat- 
tering amplitude. There is an orientation dependence only for dipoles of the size 
comparable with the impact parameter b. 
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This suggests the following picture of the "target" described by the amplitude 
eq. (|66p . It is a position dependent distribution with a center at x = 0. The dis- 
tribution itself is locally isotropic. This means that it is not anything like a single 
electric field configuration, which has locally a direction in space, but rather a locally 
isotropic ensemble of such field configurations. The mean density of this ensemble 
decreases as 1/b 2 away from the center. A small size dipole scattering on periphery 
does not care about orientation, since it again sees a locally isotropic distribution. 
A dipole whose size is of the same order as a the impact parameter probes different 
densities depending on its orientation with respect to the "bulge" of the distribu- 
tion. Thus the angular dependence here comes due to the fact that the average 
density is impact parameter dependent and it is the differences in the density that 
are probed by the dipole. 

It is of course not necessary that the function P(p) has one dominating max- 
imum. One can imagine several broad maxima, corresponding to several centers 
of density distribution. On such a target again the angular dependence appears 
through the dipoles whose size is of the order of the distance scale on which the 
density varies. The amplitudes of small dipoles are not correlated in this approx- 
imation. In other words, this mechanism of the angular correlations looks very 
different from the one we have discussed above. The difference is not just quali- 
tative. The length scale that is relevant for the mechanism of Ref.(pj6]) is not the 
correlation length of the color electric fields, but rather the correlation length of 
the density. In other words the correlations should be maximal not at k ~ Q s , but 
rather at k ~ Iq 1 , where Iq s is the correlation length of the saturation momen- 
tum in the transverse plane. Although naively one might think that the two are 
of the same order, this does not have to be case. Indeed in a model considered in 
Ref. (jnnj) the latter correlation length is parametrically larger, in the high energy 
limit Iq s ~ Qj 1 e cln ( 1 /° !s J . This mechanism would then lead to correlations at very 
small transverse momentum. 

7. Conclusions 

We have explained that rapidity and angular correlations are a very general feature 
of particle production at high energy. They are an almost automatic consequence 
of the boost invariance of the projectile wave function, provided the particle den- 
sity/scattering probability is large enough so that there is non- negligible probability 
to produce more than one particle at fixed impact parameter within the saturation 
radius 1/Q S . 

There are good reasons to expect that the factorization of both projectile and 
target averages is broken at leading order in 1 /N c in the kinematical domain relevant 
to the correlated production of particles. To study this question one certainly has 
to go beyond simple rotationally invariant solutions of the BK equation, but more 
importantly, include Pomeron loops in high energy evolution. While technically 
challenging, it would be very interesting to understand and quantify this effect. 
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We expect the correlated production to be maximal at momenta close to the 
saturation momentum Q s . It probes the hadronic wave function at transverse res- 
olutions below 1/Qs, and should be dominated by r ~ for trivial phase space 
reasons. To be able to describe fluctuations at high energy quantitatively one needs 
to understand the evolution of the high momentum modes of the target fields in the 
dilute regime and the semihard momenta in the regime of the onset of saturation. 
In a sense one has to include the KLWMIJ evolution for the high momentum modes 
and JIMWLK evolution for the low momentum ones within the same framework. 
The situation is even more complicated, since one needs to understand the structure 
of the wave function around the saturation momentum, where the Pomeron merging 
and splitting (JIMWLK and KLWMIJ contributions) are of equal importance and 
one cannot just roughly switch from one description to the other. In other words, 
the interesting momentum range is sensitive at the leading order to the Pomeron 
loop contributions. 

Current numerical implementations Ref. (fT4l I7T| l72|) do not include effects of 
KLWMIJ evolution at high momenta, nor Pomeron loops at Q Sl and as a result 
only include 1 /N c suppressed contributions to the correlated emission. In our view 
it is crucial to understand the effects mentioned above and include them in numer- 
ical calculations in order to have a theoretically robust and reliable description of 
correlated particle emission at high energy. 

Acknowledgments 

The work of AK is supported by DOE grant DE-FG02-92ER40716. The work of ML 
is partially supported by the Marie Curie Grant PIRG-GA-2009-256313 and THE 
ISRAEL SCIENCE FOUNDATION. M.L. thanks UCONN's Physics Department 
for the hospitality during which this work was initiated. 

Appendix A. Derivation of single and double inclusive production 



In this Appendix, we sketch the derivation of Eqs. (|7ll0p . We start from the wave- 
function ([2]) and follow Refs. (|62l [63]) for the derivation. Two observables of the 
interest arc 



At high collision energy, each projectile's gluon scatters off the target with the 
eikonal S'-matrix 



probabilities 




(A.l) 




(x,r,)\S\a a (x,r,)) = S ab (x) 



V exp{iT c a c t (x)} 



(A.2) 
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In the light cone gauge (A~ = 0), at(x) = J dx~ A + (x) is a field of the target on 
which the gluons scatter. To compute one of the observables O in the final state, we 
have to include initial and final state radiation as encoded in fl (see [2]), to account 
for decoherence as a result of collision process, and finally to average over all possible 
configuration of projectile and target. The result have the following form 

(6) P>T = (V>l^ f (1 - S*)Sl 6 & (1 - S)S}\ip) PtT (A.3) 

and after trivial algebra over the soft gluons we arrive for single inclusive production 
to Eq. ([7]). For the double inclusive the final result reads 

with the amplitude 

A${k,p) = f e ikz+ipu f (A.5) 

J U,Z J Xi,X2 

x{fi(z - Xl ) [S( Xl ) - S(z)] P {x 1 )} a [f j {u - x 2 ) \S{u) - S{x 2 )] p{x 2 )) h 
-f J x fi{z - n)/j(u - x x ){ [S(xi) - S(z)] p(xi) + S^xi)] Y" 

+9 I h{z-u)f 3 {u-x 1 ){{S(z)-S(u))p{x 1 )S\u)} ab . 

and we have defined p = —iT a p a . We note that in this expression the gluon with 
momentum p is assumed to have larger rapidity, and thus the emission of the two 
gluons is not completely symmetric. 

In Eq. (|10[) we keep only the term with four p. This term is leading in the limit 
of large color density p ~ l/g. One should keep in mind, h owever that in this limit 
other terms not included in Eq. (|A.6|) are equally importantP^JJ^lMI The second 
term in Eq. (|A.6j) corresponds to production of two gluons emitted from the same 
color source in the incoming projectile wave-function. The third term corresponds 
to the process whereby the softer gluon has been emitted in the wave function by 
the harder one, with both gluons subsequently produced in the collision. In terms 
of BFKL ladders, the (square of the) first term is a part of the diagram containing 
two independent ladders, while the (square of the) last two terms describe emission 
of two gluons contained in the same BFKL ladder. 
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